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1.0  INTRODUCTION 


Contamination  from  the  exhaust  plume  of  spacecraft  thrusters  can  pose 
severe  problems  for  long-life  satellite  systems.  Optical  components  (particu- 
larly if  cryogenically  cooled),  thermal  control  surfaces,  and  solar  cells  are  all 
highly  susceptible  to  reduced  performance  if  coated  with  contaminants  from 
rocket  engine  exhausts.  Analytical  tools  are  needed  to  assess  contamination 
problems  both  in  the  assessment  of  potential  contamination  problems  during  the 
spacecraft  design  phase  and  in  the  analysis  of  on-orbit  contamination  problems. 
The  CONTAM  computer  code  was  developed  by  McDonnell  Douglas  Astronautics 
Corporation  (MDAC)  under  sponsorship  of  the  Air  Force  Rocket  Propulsion 
Laboratory  (AFRPL)^^^  to  provide  these  analytical  tools  for  liquid  bipropellant 
and  monopropellant  engines. 

Although  the  CONTAM  code  has  been  in  existance  for  a number  of  years,  no 
experimental  verification  of  the  code  has  been  accomplished.  This  is  primarily 
due  to  the  lack  of  applicable  experimental  data.  Recently,  however,  the  AFRPL 
sponsored  a program  at  the  Jet  Propulsion  Laboratory  (JPL)  designed  to  obtain 
mass  deposition  rates  from  the  plume  of  a 0.  1 Ibf  monopropellant  hydrazine 
thruster.  This  data  was  then  used  by  AFRPL  in  an  attempt  to  assess  the  usabil- 
ity, accuracy,  and  deficiencies  of  the  monopropellant  portion  of  the  CONTAM 
code. 

This  report  presents  a description  of  the  AFRPL' s CONTAM  code  evaluation 
effort.  In  addition,  general  provisions  for  assessing  deposition  from  monopro- 
pellant hydrazine  thrusters  of  all  thrust  classes  are  discussed. 


(1)  R.J.  Hoffman,  ct,  aJ.  , Plume  Contamination  Effects  Prediction.  The 

CONTAM  Computer  Program,  Version  II,  AFRPL-TR-73-46,  August  1973. 
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2.  0 DESCRIPTION  OF  THE  JPL  EXPERIMENT 


2.  1 Overview 

The  data  base  for  the  evaluation  of  the  CONTAM  code  is  taken  from  mass 
deposition  rate  measurements  made  on  the  plume  of  a 0.  1 Ibf  hydrazine  thruster. 
Complete  details  of  the  JPL  experiment  are  documented  in  Reference  2.  For 
the  sake  of  completeness,  the  experiment  is  summarized  in  this  section. 

All  of  the  JPL  testa  were  conducted  in  the  JPL  MOLSINK  facility  using  a 
hydrazine  thruster  supplied  by  Hamilton  Standard.  Basically,  the  JPL  effort 
consisted  of  three  phases.  First,  deposition  rate  measurements  were  made  with 
quartz  crystal  microbalancea  (QCM's),  Next,  a series  of  about  100,000  pulses 
was  put  on  the  engine  using  a combination  of  duty  cycles.  Finally,  additional 
deposition  rate  measurements  were  made  to  assess  the  aging  effect  of  the  100,000 
pulses  on  contamination. 

The  deposition  rate  measurements  were  made  using  five  QCM's  mounted  in 
the  far  field  of  the  plume  at  angles  of  approximately  0°,  +^15°,  and  +30°  relative 
to  the  plume  centerline.  The  thruster  was  fired  with  monopropellant  grade 
hydrazine  [MIL-P-26536C,  Table  I is  a summary  of  the  JPL  deposition 

rate  tests.  Table  II  is  a description  of  the  aging  test  series  performed  by  JPL 
on  the  engine. 

2.  2 Test  Data 

The  QCM  measures  mass  deposition  by  equating  a change  in  frequency  to  a 
change  in  mass.  By  calculating  the  slope  of  the  mass  change  with  time,  one  can 
obtain  a mass  deposition  rate.  JPL's  frequency  measurements  were  made  every 
120  seconds  and  thus  the  mass  deposition  rate  measurements  are  actually  average 
rates  over  a period  of  about  twelve  pulses,  including  on  time  and  off  time. 


(2)  R,  Passamaneck  and  J,  E.  Chirivella,  "Contamination  Measurements  for  a 
0.  1 Ibf.  Monopropellant  Thruster,"  JANNAF  Ninth  Plume  Technology 
Convention,  February  1976. 
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TABLE  I 


SUMMARY  OF  JPL  DEPOSITION  RATE  TESTS 


Duty  Cycle  Bed  Temperature  QCM 

Nomenclature  sec  on/sec  off  °K  Temperature  °K 


Baseline 

0.  1/10 

478 

144, 

172, 

200 

Baseline  (with  Heat  Shield 

0,  1/10 

478 

106, 

144, 

172, 

on  Engine) 

200 

Short  Pulse 

0. 025/5 

478 

144, 

172, 

200 

Long  Pulse 

0. 2/20 

478 

144, 

172, 

200 

294°  K Cat  Bed 

0,  1/10 

294 

172 

367°  K Cat  Bed 

0. I/IO 

367 

172 

Baseline  (Water  added)* 

0.  I/IO 

478 

144, 

172, 

200 

Baseline  (aged)** 

0.  I/IO 

478 

144, 

172, 

200 

589°  Cat  Bed*** 

0.  1/10 

589 

172 

H2O  mass  fraction  = 0.  0181,  all  other  H2O  mass  fractions  = 0.  0071 

>lol<  After  100,000  pulses  (see  Table  II);  since  model  does  not  account  for  aging, 
this  case  was  not  run  analytically. 

This  case  was  not  run  analytically. 

During  the  mass  deposition  teats,  JPL  measured  thruster  wall  temperature 
to  determine  wlien  equilibrium  had  been  reached.  During  the  aging  tests,  chamber 
pressure  (Pc)  was  also  monitored.  Photographs  of  the  Pc  traces  displayed  on  an 
oscilloscope  provided  the  peak  pressures  for  the  aging  tests. 


TABLE  II 

JPL  AGING  TESTS 


; SERIES  # 

Time  on  (sec) 

Time  off  (sec) 

No.  Pulses 

1 

0.  020 

120 

50 

’ 2 

0.  035 

120 

50 

' 3 

0.  050 

120 

50 

; j 4 

0.  100 

120 

50 

i 

^■1  ' 

0.  250 

300 

5 

! 6 

0.  015 

15 

100 

j ^ 

0,  500 

120 

10 

I 8 

1, 0 

120 

10 

9 

0.  125 

12.  375 

1000 

' 10 

0.  04 

1.  16 

5000 

11 

0.  04 

0.  40 

7000 

1 

12 

1.0 

10.  0 

100 

1 

0.  10 

0.  90 

6000 

14 

0.  075 

0.  525 

6000 

' 15 

0.  090 

0,  100 

7000 

i 

0.  020 

0.  020 

1000 

( 

TOTAL’S* 

33425 

>1'  Series  1-16  is 

repeated  three  times 

for  =100,  000  total  pulses. 

3.  0 CONTAM  MODEL  IMPLEMENTATION 


3.  I Model  Deacription 

The  CONTAM  code  fur  monopropellant  hydrazine  thrusters  actually  consists 
of  a series  of  four  separate  routines.  The  N2H4  routine  calculates  the  end- 
of-bed  properties  suitable  for  input  into  the  flowfield  parts  of  the  code.  The 
MULT  RAN  routine  performs  a method-of-characteristics  calculation  of  the 
noazle  and  plume  flowfield  characteristics  following  a transonic  analysis  through 
the  throat.  A one-dimensional  chemical  kinetics  calculation  along  the  MULTRAN, 
pressure-defined  streamlines  is  then  performed  with  the  KINCON  routine.  This 
routine  also  includes  the  capability  to  predict  plume  condensation  of  a single 
species.  The  resultant  velocities,  densities,  and  species  concentrations  along 
streamlines  are  used  to  calculate  the  mass  flux  points  within  the  plume.  The 
SURFACE  routine  calculates  deposition  ra  ■'s  to  the  spacecraft  surface  and 
resultant  changes  in  surface  properties.  Since  the  emphasis  of  this  work  was 
placed  on  mass  deposition  rates  as  a function  of  surface  temperature,  this 
routine  was  not  exercised.  When  the  CONTAM  code  is  used  for  bipropellant 
engines,  a transient  combustion  chamber  code  (TCCl  ip  substituted  for  N2H4. 
while  the  other  codes  remain  the  same.  The  data  base  for  this  effort  was  a mono- 
propellant engine  and  thus  TCC  is  not  discussed  in  this  report, 

Since  SURFACE  was  not  used,  the  net  deposition  rate  predicted  at  the  QCM 
surface  had  to  be  calculated  by  subtracting  the  evaporation  rate  of  each  species 
from  the  plume  deposition  predicted  by  N2H4,  MULTRAN.  and  KINCON.  Fol- 
lowing is  a discussion  of  each  model  and  of  the  technique  used  in  calculating  net 
deposition  rate. 

3.  2 N2H4  Routine 

The  N2H4  code  is  a one-dimensional,  non-steady  state  model  which  calcu- 
lates end-of-bed  properties  using  finite  difference  solutions  of  mass  transfer, 
heat  transfer  and  chemical  reaction  rate  equations.  In  actuality,  N2H4  is  taken 
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with  only  minor  modification  from  the  transient  model  developed  for  NASA  by 
A,  S,  Kesten,  A very  complete  discussion  of  the  equations  and  limi- 

tations used  in  the  transient  model  is  detailed  in  the  references.  The  reports 
also  include  the  derivation  or  references  for  the  physical  properties  used  in  the 
model. 

Table  III  is  a listing  of  the  input  parameters  used  in  the  "baseline"  test  case. 
The  steady-state  parameters  are  obtained  from  a steady-state  model  developed 
by  Kesten.  The  geometrical  inputs  for  the  steady  state  model  are  the 

same  as  those  for  N 21-14.  Since  N2H4  estimates  by  extrapolation  to  steady- state 
conditions,  the  running  of  the  steady  state  model  is  essential.  It  is  not,  however, 
included  in  the  CONTAM  code  and  must  be  obtained  separately. 

During  the  course  of  running  N2H4  several  difficulties  were  encountered. 

The  most  serious  of  these  was  the  way  in  which  N2H4  modeled  the  thruster's 
thermal  characteristics.  I'r2H4  does  not  contain  any  provisions  for  a heat  shield 
nor  for  cooling  to  a space  temperature  different  than  the  temperature  to  which 


(3)  A.  S.  Kesten,  Analytical  Study  of  Catalytic  Reactors  for  Hydrazine 
Decomposition.  United  Aircraft  Research  Laboratories  Report 
G910461-24,  Second  Annual  Progress  Report,  Contract  NAS-7-458,  May 
1968. 

(4)  A.  S.  Kesten,  Analytical  Study  of  Catalytic  Reactors  for  Hydrazine 
Decomposition.  United  Aircraft  Research  Laboratories  Report 
H910461-38,  Third  Annual  Progress  Report,  Contract  NAS-7-458,  May  1969. 

(5)  D.  B.  Smith,  E,  J.  Smith,  and  A.  S.  Kesten,  Analytical  Study  of  Catalytic 
Reactors  for  Hydrazine  Decomposition.  United  Aircraft  Research 
Laboratories  Report  H910461-37,  Computer  Program  Manual,  Transient 
Model,  May  I969. 

(6)  A,  S.  Kesten,  Analytical  Study  of  Catalytic  Reactors  for  Hydrazine 
Decomposition.  United  Aircraft  Research  Laboratories  Report  F91 046 1-12, 
First  Annual  Progress  Report,  Contract  NAS-7-458,  May  1967. 

(7)  D.  B.  Smith,  E,  J.  Smith,  and  A.  S.  Kesten,  Analytical  Study  of  Catalytic 
Reactors  for  Hydrazine  Decomposition.  United  Aircraft  Research  Labora- 
tories Report  G910461 -30,  Computer  Programs  Manual,  One-Dimensional 
and  Two-Dimensional  Steady  State  Models,  August  1968, 
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TABLE  III 


N2H4  INPUTS.  BASELINE  CASE 


QUANTITY 

DESCRIPTION 

VALUE 

UNITS-'- 

Su-;uly 

-Slate  Parameters 

GOSS 

Steady  State  Mass  Flow  Rate 

0. 1 1772 

1 bm  / ft  ^ - s ec 

MRSS 

Steady  State  Molecular  Weight 

12.  525 

Ibm/lb-molc 

FSS 

Steady  Slate  Chamber  Pressure 

81. 04 

psia 

TSS 

Steady  State  Temperature 

2037 

“r 

TVAP 

Hydrazine  Vaporization  Temp. 

804 

“r 

llIV 

Enthalpy  of  Liquid  Vapor 

709.  2 

BTU/lbm 

Operational  Variables 

F' 

Interface 

f 

Flow  Rate  Through  Buried  Inj. 

98.  1 

Ibm/ft^-si-c 

HF 

Enthalpy  of  Feed 

78. 

BTU/lbm 

PF 

Pressure  of  Feed 

114. 

psia 

TA 

Ambient  Temperature 

860. 

"r 

TPl 

Initial  Catalyst  Temp«-rature 

860. 

"r 

TLIQ 

Initial  Concentrations  of 

N2II4,  NH3,  H2,  N2 
Tentperaltirc  of  Liquid 

1.  OxlO"'-* 
860. 

Ibm/ft^ 

"R 

V ' 

TS 

Hydrazine 

Background  Temperature 

36 

‘’r 

PCI 

Duty  Cycle 

Initial  Chamber  Pressure 

. 1 sec  on/ 

1 0 S('c  off 
0.  001 

psia 

Clhambt 

■r  t;  Catalyst  Characteristics 

AC  Cross  Se<'tional  Area  of 

3.  98x10"* 

ft*^ 

MW 

React  ion 

Chamber  Wall  Thermal  Mast 

0. 01378 

Ibm 

VC 

Free  Volume  of  Chamber 

4.  95x10'*' 

ft  3 

AW 

Chambi-r  Wall  Specific  Heat 

0.  1 

BTU/lbm-''R 

HA 

F'orced  Convection  Heal 

0. 

BTU/fl^- 

IlAl 

Tr.ansfer  Coefficient 

Natural  Convention  Heat 

0. 

sec-“R 

BTU/fl^- 

1IA2 

T ran s for 

Radiative  Heat  Transfer 

1. 5x10-' 3 

sec- 1 . 25-"! 
BTU/fl^- 

AS 

Coefficient 

I leat  Shield  Area 

(No  shieUl) 

sec-^R"* 
ft  2 

/, 

Length  of  Bed 

0. 0358 

ft 

/.0 

Axial  Distance  to  Enel  of 

0. 0108 

ft 

KADI' 

Buried  Injector 

Radius  of  Catalyst  Particles 

8.  '>8x10"* 

ft 

SUK  F T 

|2S  30  Mesh) 

Catalyst  Particle  Surfae c 

2070 

ft  “^  / ft  3 

VOID'!' 

Area/Unit  Volume  of  Bed 

1 111  e rpa  rt  ic  li-  Void  Fraction 

0.  38 

Qii.inl  It  ii-s  arc  in  iinilH  required  for  propram  inpol. 
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the  catalyst  bed  is  heated.  These  provisions  were  added  by  AFRPL.  since  they 
were  required  for  modeling  the  JPL  tests.  In  addition,  original  estimates  for 
natural  convection  and  radiation  coefficients  gave  a poor  thermal  fit  to  the  base- 
line data.  Since  the  purpose  of  this  effort  was  to  predict  deposition  and  not  wall 
temperature,  it  was  decided  to  approximate  the  baseline  thermal  response 
empirically  using  only  radiation.  A radiation  coefficient  which  gave  a close  fit 
to  the  baseline  thermal  response  was  found  (Table  III)  and  used  throughout. 

Figure  1 shows  the  baseline  fit. 

When  N2H4  is  operated  in  the  pulse  mode,  the  peak  pressure  predicted  for 
a given  pulse  is  higher  than  that  for  the  previous  pulse.  This  difference  decreases 
gradually  until  finally  each  pulse  is  predicted  to  be  identical  to  the  next,  and  for 
the  purpose  of  comparison  with  experiment,  is  considered  to  be  "steady- state". 
Figure  2 shows  a comparison  of  the  pressure  traces  for  pulses  1,  5,  45,  and  60 
of  the  baseline  case.  After  pulse  45,  the  traces  are  nearly  identical.  Upon 
re-examining  Figure  1,  it  can  be  observed  that  the  temperature  also  levels  off 
after  about  45  pulses  (450  seconds).  Past  pulse  60,  the  model  predicts  no  change 
in  engine  parameters,  and  it  is  these  steady-state  parameters  that  were  used  to 
develop  input  for  the  subsequent  contamination  prediction  routine.  However, 

N2H4  predicts  values  as  a function  of  time  throughout  a pulse  while  MULTRAN 
and  KINCON  require  properties  at  a discrete  time.  The  simplest  technique  is 
to  average  the  parameters  over  the  entire  pulse.  For  chamber  pressure  and 
temperature  the  average  is  obviously  an  acceptable  approach.  In  the  "steady- 
state"  regime  (past  pulse  60)  the  temperature  is  fairly  constant.  Also,  the 
pressure  rise  and  fall  are  so  fast  that  the  chamber  pressure  is  essentially  at 
the  average  for  the  entire  pulse.  Conversely,  the  end-of-bed  species  concentra- 
tions, essential  for  calculation  of  mass  flux,  vary  quite  a bit  and  thus  the 
accuracy  of  the  averaging  technique  ior  these  properties  has  to  be  evaluated. 

The  validity  of  this  technique  is  discussed  in  the  Transient  Effects  section  (3.  6). 
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Figure  1.  Baseline  Thermal  Profile  Matching 


Table  IV  shows  the  inputs  to  MULT  RAN  and  KINCON  that  were  calculated 


by  N2H4.  Since  the  number  of  pulses  for  the  model  to  reach  steady  state  varied 
with  conditions,  the  pulse  used  for  calculations  is  also  indicated.  It  should  be 
noted  that  the  "water  added"  inputs  were  taken  from  the  baseline  case  and 
adjusted  for  the  higher  mass  fraction  of  water  in  the  propellant. 


TABLE  IV 

N2H4  - CALCULATED  INPUTS  TO  MULTRAN  AND  KlNCON'^t 


CASE  (Pulse 
for  steady 
state] 

Specific 

Heat 

Ratio 

Chamber 

Pressure 

PSIA 

Chamber 

MASS  FRACTION* 

R 

NH3 

Nz 

Hz 

N2H4 

Baseline  (6O) 

1. 1941 

59.  46 

1225 

. 51354 

. 46365 

. 01571 

4x10*^ 

Short  Pulse  (60) 

1.  1876 

22.  47 

1009 

. 52128 

, 46023 

. 01 1.39 

1.  7x1 0 

Long  Pulne  (40) 

1.2131 

66.  80 

1225 

.45484 

. 51691 

.02115 

1. 0x10'^ 

70°  Bed  (40) 

1.  1921 

58.  96 

1143 

. 53.182 

. 44424 

. 01684 

2.  3x10“^ 

200°  Bed  (60) 

1. 1917 

59.  17 

1178 

. 53003 

. 44673 

.01614 

2.  5x10''^ 

Water  Added 

Same  as 
increase 

Baseline  e 
water 

xcept  H 

2O  Fluj 

was  ac 

justed  i 

or  1% 

+ quantities  were  averaged  over  tota.l  pulse 

4=  quantities  are  in  units  required  for  program  input 

mass  fraction  of  H^O  = . 0071  except  for  "water  added"  case,  where  H2O  mass 
fraction  = .0181 


3.3  MULTRAN  Routine 

The  MULTRAN  routine  of  CONTAM  performs  a method  of  Characteristics 
(MOC)  calculation  following  a transonic  analysis  through  the  tliroat,  Inputs  to 
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the  MULTRAN  code  from  N2H4  have  already  been  discussed  in  Section  3.  2,  The 
engine  design  parameters  used  for  MULTRAN  input  are  presented  in  Table  V. 
The  actual  radius  of  curvature,  2.  0,  had  to  be  changed  to  2.  667  since  the  code 
was  not  designed  to  handle  small  radii  of  curvature.  From  a sensitivity  study 
of  this  parameter,  varying  it  over  the  range  that  the  code  would  allow,  it  can  be 
concluded  that  the  results  are  insensitive  to  this  change. 


TABLE  V 

MULTRAN  INPUTS 


Quantity 

Description 

Values 

Units+ 

THID 

Inlet  Angle 

30 

Degrees 

RT 

Throat  Radius 

. 00125 

Feet 

RRT>:< 

Throat  Radius  of 
Curvature 

2.  667 

Normalized  to  Throat 
Radius 

EPS 

Nozzle  Area  Ratio 
(Core) 

55.  3 

- - 

* Actual  Value 

- 2.  0,  see  text. 

+ Quantities  are  in  units  required  for  program  input. 

In  the  operation  of  MULTRAN.  it  was  found  that  the  model  would  not  calcu- 
late values  beyond  a distance  of  30  - SO  throat  radii,  nor  would  the  model  accept 
a Prandtl-Meyar  expansion  angle  of  greater  than  70°.  Upon  completion  of  the 
MOC  calculation,  the  code  sets  up  a series  of  pressure  defined  streamlines  that 
are  subsequently  read  into  KINCON.  Figure  3 shows  a typical  plot  of  MULTRAN 
streamlines.  If  the  99.99%  streamline  of  that  figure  is  examined,  it  can  be  seen 
that  the  points  (shown  by  clrclesjdo  not  all  fall  on  the  expected  line.  Thus  the  true 
angle  that  the  streamline  would  strike  a surface  a large  distance  away  from  the 
thruster  could  not  be  accurately  determined.  This  plus  the  inability  of  the  model 
to  expand  the  plume  beyond  70°  had  a great  effect  on  the  off-centerline  calculations. 
This  effect  will  be  discussed  further  in  Sections  4 and  5. 
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MULTRAN  Generated  Streamlines 


3.  4 KINCQN  Routine 


KINCON  was  used  as  described  in  the  CONTAM  manual,  with  the 
pressure-defined  streamlines  from  MULT  RAN  and  species  mass  fraction  pro- 
duced by  N2H4  as  input,  with  the  following  exceptions:  (a)  hydrazine  was  not 
included  in  the  KINCON  calculation  because  data  for  it  is  not  included  in  the 
standard  KINCON  thermodynamic  library;  (b)  water  was  added  to  the  flow  field 
in  the  same  amount  as  was  found  in  the  propellant.  (The  nitrogen  concentration 
was  adjusted  slightly  to  make  the  mass  fractions  sum  to  one.  ) No  chemical 
reactions  were  included  in  the  calculation,  except  for  the  dummy  reaction 
NHj— ^■NHj  (the  program  requires  at  least  one  reaction).  The  species  mass 
fractions  were  thus  frozen  at  the  chamber  exit  values.  This  approximation 
should  be  quite  good  because  of  the  relatively  low  pressures  and  temperatures 
in  the  nozzle  and  plume  compared  with  the  chamber  conditions.  Eventhough 
hydrazine  was  not  included  as  part  of  the  flow  field  calculation  itself,  it  was 
added  back  into  the  composition  at  its  end-of-bed  value  for  the  mass  deposition 
rate  calculation.  Although  KINCON  also  has  the  capability  of  calculating  con- 
densation of  one  species  in  the  flow  field,  this  option  was  not  exercised  for  this 
work. 

3.  5 Calculation  of  Net  Deposition  Rates 

The  output  from  KINCON  was  used  to  obtain  mass  fluxes  at  points  along  each 
streamline  from  the  equation: 

rh(r)  = p(r)  • V(r)  (I) 

where:  p(r)  = the  gas  density  at  r in  g/cm^ 

V(r)  = the  gas  velocity  at  r in  cm/sec 

r = the  distance  along  the  streamline  from  the  throat  in  cm 
rh(r)  a the  mass  flux  per  unit  area  at  r in  g/cm  - sec 
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This  mass  flux  per  unit  area  can  be  converted  to  a mass  flux  per  vinit  solid 
angle  via  the  equation: 


drh 

dlT 


(r)  s m(r) 


’cos  </> 


(2) 


where: 


= the  mass  flux  per  unit  solid  angle  in  g/sec-steradian 


= the  angle  between  the  direction  of  flow  and  the  normal  of  the  surface 
through  which  the  flux  is  calculated. 

In  principle,  rh(r)  can  be  calculated  at  positions  along  any  streamline  to  a 
distance  where  the  QCM's  are  located  to  compare  with  the  QCM  data.  In 
practice,  It  is  found  that  becomes  constant  at  relatively  close  distances 

from  the  exit  plane  compared  with  the  distance  to  the  QCM's,  thus  demonstrating 
that  the  plume  far  field  is  similar  to  a source  flow.  The  calculation  with 
KINCON  is  only  done  out  to  the  point  where  becomes  constant.  The  mass 

flux,  rh(r),  can  be  found  by  solving  equation  (2)  at  any  location  r.  The  angle^  was 
taken  to  be  0“  since  the  QCM  data  had  been  corrected  for  the  inclination  effect. 
Once  m(r)  at  the  QCM  is  calculated,  the  deposition  rate  for  each  contaminant 
species  is  calculated  from: 


rh^  = X^m 


(3) 


where:  rh^  = the  mass  deposition  rate  for  species  i in  g/cm^  - sec 
= the  mass  fraction  of  species  i in  the  flow  field 
m = the  total  mass  flux  per  unit  in  g/cm^  - sec 


The  mass  flux  predicted  by  CONTAM  is  the  flux  for  the  total  run  time;  the 
time  for  which  the  chamber  pressure  was  above  its  initial  value.  To  calculate 
the  mass  flux  rate  hitting  the  crystal,  it  is  necessary  to  average  the  flux  over 
the  sum  of  the  on  and  off  time  (See  Appendix  A), 
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The  evaporation  rate  for  each  species 
from  reference  8; 


= (5.83  X 10-2)  YPy 


. was  calculated  using  the  expression 


(4) 


where;  rhg,  = the  evaporation  rate  in  g/cm^  - sec  of  species  i 
Y is  an  accomodation  coefficient  (assumed  = 1) 

Pv  = equilibrium  vapor  pressure  in  Torr  of  species  i at  temperature  T 
MW  = molecular  weight  of  species  i in  g/mole 
T = QCM  temperature  in  °K 


The  choice  of  unity  as  the  value  of  the  accomodation  coefficient  is  equivalent 
to  assuming  that  any  deposit  equilibrates  to  the  QCM  temperature.  Also,  this 
expression  assumes  zero  external  pressure  and  is  not  valid  for  external  pres- 
sures of  about  the  same  magnitude  as  the  vapor  pressure  or  larger.  If  this 
condition  occurs,  the  actual  evaporation  rate  will  be  smaller  than  that  calculated 
with  equation  4.  However,  in  the  experiment  where  the  vacuum  chamber  was 
constantly  cryopumped,  equation  4 should  be  valid. 

IQ) 

From  the  CRC  Handbook  of  Chemistry  and  Physics^  , the  vapor  pressure 
for  solid  ammonia  was  represented  by: 

_ -1630. 1 5 _ ,»> 

log  + 9.9974  (5) 


where:  Py  = ammonia  vapor  pressure  in  torr 
T = Temperature  in  °K 


(8)  R.J.  Hoffman,  et.  ^1.,  Op.  Cit.  , page  230. 

(9)  Handbook  of  Chemistry  and  Physics,  The  Chemical  Rubber  Company, 
Cleveland,  Ohio,  19^7  - 68  Edition. 


18 


This  equation  is  an  empirical  curve  fit  valid  for  temperatures  between  146°K 
and  195'^K.  Ice  vapor  pressure  data  was  obtained  from  the  equation: 


Log  Py  = + 3.  5721988  log  T -3.  97203  x 10'^  x T -1.  7649  x 10'^ 


t2  + 1. 901973 


(6) 


where:  = ice  vapor  pressure  in  torr 

T a temperature  in  °K 

This  equation  was  obtained  from  a theoretical  derivation^  which  has  been 

shown  to  match  experimental  data  between  173°K  and  373°K  to  within  the 

experimental  accuracy.  Hydrazine  vapor  pressure  was  obtained  from  a cor- 

B 

relation  of  vapor  pressure  using  a three  constant  expression,  In  P^  = A" 

XtC 

where  C is  an  empirical  function  of  the  critical  temperature  and  A and  B are 
calculated  from  known  vapor  pressures.  The  equation  is  as  follows: 

3273 


‘"Pv=  ^5 


(7) 


The  net  deposition  for  each  species  is  then  calculated  as 


mnet  i = nii  - mg. 


(8) 


If  linnet  i then  rhjj^^  ^ = 0,  The  total  net  deposition  is  calculated  by  summing 
over  all  species  considered. 


N 

’^net  ~ E *^net  i 
i=l 


(9) 


For  the  purpose  of  this  analysis  only  N2H4,  NH3,  and  H^O  are  considered. 
The  evaporation  rates  of  K2  and  N2,  the  other  species  of  N2H4  decomposition  are 
too  high  to  deposit  at  the  crystal  temperatures  used  in  this  analysis. 


(10)  G.  JancBO,  et.  , J.  Phys.  Chem.  , 1970,  74(15),  2984-9. 
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3.  6 Transient  Effects 


As  mentioned  in  the  proceeding  sections,  the  N2H4  model  generates  inputs 
as  a function  of  time  during  a pulse.  To  use  MULT  RAN  and  KIN  CON,  it  was 
necessary  to  average  these  inputs  over  the  pulse  duration.  The  validity  of  this 
technique  was  assessed  by  comparing  the  average  of  mass  fluxes  calculated  for 
several  times  during  a pulse  with  the  mass  flux  calculated  from  average  pulse 
properties.  Based  on  the  results  of  this  comparison,  (see  Figure  4),  only  one 
set  of  calculations  was  needed  for  each  pulse.  A similar  check  of  the  individual 
mass  fluxes  for  each  specie  also  indicated  that  the  pulse  average  quantity  is  a 
reasonable  value,  except  perhaps  for  hydrazine  which  has  the  most  fluctuation 
during  a given  pulse.  The  approximation  in  all  cases  can  be  expected  to  be 
better  for  longer  pulses  when  the  average  is  closer  to  the  statistical  mode. 

As  the  change  from  pulse  to  pulse  of  N2H4  predicted  parameters  begins  to 
become  minimal  (see  Section  3,  2),  so  too  do  the  differences  in  CONTAM  predicted 
mass  deposition  rates.  Figure  5 indicates  this  for  the  baseline  case,  It  can  be 
seen  that  there  is  a gradual  decline  in  the  deposition  rate  over  the  first  10-15 
pulses.  By  pulse  60,  the  deposition  rate  has  leveled,  and  a "steady  state"  is 
obtained.  This  deposition  rate  was  used  to  calculate  the  net  deposition  rates 
for  comparison  with  JPL  data.  Since  JPL  measured  QCM  frequency  only  every 
120  seconds  (approximately  12  pulses),  it  was  impossible  to  compare  the 
CONTAM  predictions  during  early  pulses  with  JPL's.  In  fact,  JPL's  data  comes 
from  readings  taken  after  at  least  100  pulses.  In  making  the  comparisons,  the 
initial  CONTAM  transient  must  be  ignored.  Since  changes  in  deposition  rate  and 
not  absolute  mass  are  being  measured,  this  limitation  should  not  be  serious. 


4.  0 PRESENTATION  OF  RESULTS 


4.  1 N2H4  Data  Comparison 

Although  the  comparison  between  JPL  measured  and  CONTAM  predicted 
mass  deposition  rates  is  of  primary  concern  in  this  effort,  the  engine  data  taken 
by  JPL  during  the  tests  (see  Section  2.  2)  offered  a chance  to  evaluate  the  output 
from  the  N2H4  code. 

Table  VI  shows  a comparison  between  experimental  and  predicted  values  of 
peak  pressure  for  those  aging  duty  cycles  (see  Table  II)  for  which  N2H4  was  run. 
It  should  be  remembered  that  JPL's  pressure  value  had  to  be  interpolated  from 
photographs  of  oscilloscope  traces.  Thus,  the  agreement  for  all  but  series  No.  1 
is  reasonable.  Based  on  the  short  on-time,  one  must  assume  that  the  reported 
experimental  value  for  series  No.  1 is  in  error. 

Table  VII  presents  the  maximum  and  minimum  temperatures  recorded  by 
JPL  for  six  of  the  aging  series  (see  Table  II)  and  the  model  predictions  for  these 
values.  The  correlation  between  the  predicted  and  experimental  maximum 
temperatures  is  good,  (0.  96)  although  for  long  off  times  (>120  seconds)  or  for 
long  on  times  (>.15  second)  the  final  wall  temperature  is  overestimated.  This 
is  understandable  if  one  considers  the  empiricism  used  to  model  wall  cooling 
(see  Section  3.2).  Figure  6 compares  several  N2H4  predicted  wall  temperature 
profiles  with  JPL  measured  data.  When  the  particular  duty  cycle  was  close  to 
the  baseline  (0.  I on/ 10  off)  the  fit  was  good  (see  Figure  1),  As  the  duty  cycle 
varied  from  the  baseline  case  the  fit  became  poorer. 

4.  2 Centerline  Mass  Deposition  Rates; 

The  centerline  mass  deposition  rates  predicted  by  CONTAM  for  all  of  the 
JPL  cases  modeled  are  presented  in  Table  VIII.  It  should  be  remembered  that 
these  are  the  predicted  rates  at  which  mass  strikes  the  QCM.  To  calculate  the 
net  deposition  rate,  one  must  subtract  the  evaporation  rate  of  each  species  from 
its  deposition  rate;  tlie  values  shown  in  Table  VIII  are  independent  of  QCM 
temperature.  Appendix  A presents  a sample  calculation. 
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TABLE  VI 


PEAK  PRESSURE  COMPARISON 


Series  # 

Duty  Cycle 
on/ off  (sec) 

Experimental  Peak 
Pressure  (PSIA) 

Predicted 

Peak  Pressure  (PSIA) 

1 

. 02/120 

95.  0 

70.  0 

2 

. 03B/120 

72.  5 

71. 2 

3 

. 05/120 

80. 

71.4 

4 

. 10/120 

72.  5 

71.  8 

5 

. 25/300 

72.  5 

73.  6 

TABLE  VII 

TEMPERATURE  RANGE  COMPARISON 


Series  # 

Duty  Cycle 

Experimental  (°R) 

Predicted  (°R) 

MAX 

MIN 

MAX 

MIN 

1 

02/120 

873 

860 

906 

861 

2 

. 035/120 

880 

860 

904 

860 

3 

. 050/120 

881 

860 

908 

860 

4 

. 10/120 

891 

860 

889 

862 

5 

. 25/300 

923 

860 

964 

900 

8 

1. 0/120 

(1)>I<  1038 

860 

1056 

999 

(2)>l<  1043 

889 

1152 

1065 

(3)>l'  1045 

890 

1203 

1094 

o 

o 

897 

1226 

107  5 

(5)>1<  1052 

899 

1236 

Pulse  number 


^4 


BASELINE  (WITH  SHIELD) 


Experimental  Vs.  Calculated  Thermal  Profiles 


TABLE  VIII 


CALCULATED  CENTERLINE  MASS  DEPOSITION  RATES 


ENGINE  CONDITION 

DEPOSITION  RATES  AT 

QCM  (g/cm^/sec) 

Total 

NH3 

H2O 

N2H4 

xio’^ 

xlO'' 

xl09 

xl0^° 

Baseline 

9.  50 

4.  87 

6.  74 

3.80 

Baseline  (w/ shield) 

9.  33 

4.  79 

6.  67 

3.  76 

Short  Pulse 

7.  98 

4.  15 

5.  66 

1.  35 

Long  Pulse 

8.  81 

4.00 

6.  23 

8.79 

70°F  Catalyst  Bed 

9.  66 

5.  13 

6.  84 

2.  22 

2O0°F  Catalyst  Bed 

9.  56 

5.  07 

6.  79 

2.  39 

Baseline  (H2O  added) 

9.49 

4.87 

17.  73 

3.  81 

From  Table  VIII  it  can  be  observed  that,  although  mass  deposition  from  the 
short  pulse  and  long  pulse  trains  was  predicted  to  be  slightly  less  than  that  of 
the  baseline  case,  total  mass  deposition  does  not  vary  a great  deal  with  engine 
conditions.  The  fluctuations  observed  in  the  hydrazine  flux  are  due  almost 
entirely  to  fluctuations  in  hydrazine  mass  fraction  (see  Table  IV).  There  is 
little  difference  in  the  deposition  rates  of  water  and  ammonia  for  the  different 
engine  conditions,  except  of  course  for  the  "water  added"  case  in  which  the  water 
mass  fraction  is  higher.  The  order  of  magnitude  differences  in  the  mass  deposi- 
tion rates  of  the  three  species  are  directly  relatable  to  the  order  of  magnitude 
differences  in  their  mass  fractions  in  the  plume. 

4.  3 Comparisons  with  JPL  Data;  (Centerline) 

Table  IX  presents  a comparison  between  the  predicted  net  deposition  rates 
and  the  deposition  rates  measured  by  JPL. 

From  this  table,  it  can  be  seen  that,  except  for  the  200°K  case,  the  calcu- 
lated net  deposition  rates  are  all,  higher  than  the  rates  measured  by  JPL.  Both 
sets  of  data  are  consistent,  however,  in  that  deposition  rate  decreases  with 
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increasing  temperature.  This  is  to  be  expected  since  the  evaporation  rate  is  an 
extremely  strong  function  of  temperature.  Figure  7 is  a plot  of  the  predicted 
deposition  vs.  temperature  for  the  various  cases  (the  baseline  with  shield  and  the 
294°  K catalyst  bed  conditions  were  not  plotted  because  of  their  similarity  with 
other  cases).  An  examination  of  Figure  7 indicates  discrete  plateau  regions 
where  the  predicted  flux  is  constant  over  a range  of  temperatures.  In  other 
regions,  the  balance  between  deposition  and  evaporation  is  such  that  a small 
change  in  temperature  results  in  a substantial  change  in  net  deposition  rate.  By 
considering  each  of  these  areas,  one  can  get  an  excellent  feel  for  what  the  pre- 
dicted results  mean.  Above  190°K  no  net  deposition  is  predicted  for  any  case. 
This  is  because  the  evaporation  rate  of  all  the  species  is  higher  at  this  tempera- 
ture than  any  of  the  deposition  rates.  JPL.  did  measure  some  deposition  rates, 
but  this  may  well  be  a species  other  than  the  three  considered.  Compounds  from 
the  carbonaceous  impurities,  especially  aniline,  are  quite  possible. 

As  the  surface  temperature  is  lowered  from  190°  to  175°K,  the  evaporation 
rate  and  deposition  rate  of  the  hydrazine  become  balanced.  The  plateau  region 
between  160°  and  17S°K  is  thus  the  temperature  where  all  of  the  hydrazine  in  the 
flow  which  hits  the  crystal  is  expected  to  stick.  The  evaporation  rates  of  water 
and  ammonia  at  these  temperatures  still  exceed  their  deposition  rates.  In 
general  the  predictions  at  172°K  are  from  two  to  eight  times  higher  than  the  data. 
If  one  considers  the  averaging  techniques  used  to  calculate  hydrazine  flux,  and 
the  small  quantities  of  hydrazine  predicted  by  the  model,  the  lack  of  good  agree- 
ment in  the  "all  hydrazine"  region  is  not  surprising. 

As  the  temperature  of  the  surface  is  lowered  below  IbO'^K,  the  water 
evaporation  rate  quickly  becomes  comparable  to  the  water  deposition  rate.  Since 
there  is  much  more  water  in  the  flow  than  there  is  hydrazine,  a sharp  increase 
in  deposition  is  observed.  At  150°K,  the  water  deposition  rate  is  much  larger 
than  thfr  evaporation  rate,  and  another  temperature  insensitive  plateau  region  is 
reached.  Between  110°K  and  145°K,  all  the  water  and  hydrazine  which  hits  the 
QCM  is  predicted  to  stick.  At  144‘^K,  the  model  predictions  are  still  higher 
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Figure  7.  Net  Deposition  Rate  Vs.  Temperature 


than  the  JPL  data,  but  only  by  a factor  of  three  or  four.  Ironically,  one  of  the 
worst  cases  for  agreement  Is  the  water  added  case. 

Near  11Z°K  the  ammonia  begins  to  deposit.  Since  ammonia  is  present  in 
the  largest  quantity  of  the  three  species,  once  again  there  is  a sharp  rise  in  the 
net  deposition  curve,  For  the  single  106°K  comparison,  the  predicted  value  was 
a factor  of  five  higher  than  the  data.  The  results  at  this  temperature  are  thus 
consistent  with  the  other  data. 

4.  4 Off-Centerline  Mass  Deposition  Rates 

Table  X shows  a comparison  between  CQNTAM  predicted  mass  fluxes  for  the 
centerline  and  the  streamline  that  contains  99*  99%  of  the  mass  between  it  and  the 
centerline.  These  mass  fluxes  have  been  divided  by  the  distance  to  the  QCM 
squared  (see  Appendix  A)  but  have  not  been  adjusted  for  angle  from  the  center- 
line.  Unfortunately,  as  indicated  in  section  3.  3,  it  was  difficult  to  ascertain  the 
angle  at  which  a given  streamline  would  impinge  upon  the  QCM  surface, 

TABLE  X 

COMPARISON  OF  MASS  FLUX  ALONG  CQNTAM  STREAMLINE 
ENGINE  CONDITION  MASS  DEPOSITION  RATES  AT  QCM  (g/cm^/sec)(x  10''^) 


Centerline 

99.  99%  Streamline 

Ratio 

Baseline 

9.  50 

19.  7 

2.  07 

Short  Pulse 

7.  98 

7.  77 

. 97 

Long  Pulse 

8.  81 

9.  07 

1. 02 

70°F  Cat  Bed 

9.  66 

8.  03 

. 83 

200°F  Cat  Bed 

9.  56 

7.  62 

- 79 

The  99.99%  mass  fluxes  can  be  converted  to  net  deposition  rate  (by  subtracting 
the  evaporation  rate  of  each  species)  and  then  adjusted  to  the  QCM  off-centerline 
angle  (9  (multiplying  by  cos^  $).  Table  XI  presents  a comparison  between  these 
adjusted  values  and  JPL'a  measurements. 
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li>  all  cases  the  99.99%  streamline  over  estimates  the  deposition  observed 
by  JPL.  In  fact,  if  this  streamline  were  considered  to  be  at  70°,  the  maximum 
Prandtl-Meyer  expansion  angle  that  could  be  used,  tlve  predicted  deposition  is 
still  higher  than  the  JPL  15°  measurements.  This  example  indicates  that 
MULT  RAN  predicts  the  mass  to  be  concentrated  in  too  small  an  angle  and  points 
out  a serious  deficiency  in  the  code. 


5,  0 DISCUSSION  OF  RESULTS 


From  the  results  of  Table  IX,  it  can  readily  be  seen  that  the  CONTAM 
model  coupled  with  simple  evaporation  rate  calculations  predicts  centerline 
mass  deposition  in  excess  of  those  measured  by  JPL.  Such  overprediction  can 
occur  as  a result  of  either  an  over  estimation  of  the  mass  flux  arriving  at  the 
surface,  an  underestimation  of  the  evaporation  rate,  or  a combination  of  both 
these  phenomena.  JPL  data  was  taken  only  once  every  twelve  pulses,  and  thus 
deposition  and  evaporation  could  not  be  separated  experimentally.  However,  in 
this  section,  these  sources  of  error  will  be  uncoupled  as  much  as  possible  in 
order  to  evaluate  the  individual  routines  and  analysis  steps  used  in  this  effort. 
Appendix  A presents  a sample  calculation  which  should  be  helpful  in  following 
the  interactions  of  the  various  model  predictions. 

5.  1 Deposition  Rate  Predictions  Errors 


Although  the  deposition  rate  predicted  by  CONTAM  is  made  up  of  two  values, 
the  mass  flux  from  the  engine  and  the  mass  fraction  of  the  individual  species, 
it  is  possible  to  separate  these  using  the  calculations  for  the  same  engine  condi- 
tions as  a function  of  QCM  temperature.  From  Figure  7 it  can  be  seen  that  there 
are  distinct  temperature  regions  in  which  one  can  predict  which  of  the  three 
species  considered  will  deposit.  The  order  of  magnitude  changes  predicted  by 
CONTAM  for  the  106,  144,  and  17Z°K  QCM  temperature  were  substantiated  by 
the  JPL  data.  Thus,  one  c.an  postulate  with  some  certainty  that  at  172°k  only 
hydrazine  will  deposit,  at  144°K  only  water  and  hydrazine,  and  at  106”K,  water, 
hydrazine,  and  amm.onia.  Furthermore,  because  of  the  differences  in  mass 
fractions  for  the  three  species,  the  deposition  at  106°K  is  principally  amnionia, 
and  at  144‘^K  principally  water.  Since  the  mass  fraction  of  water  in  tho  pro- 
pellant is  known  with  reasonable  accuracy,  one  can  use  the  144'^K  crystal  tempera- 
ture to  evaluate  the  mass  flux  predicted.  At  this  temperature,  CONTAM 
overpredicts  the  JPL  data  by  1.  5 to  almost  5.0,  Similarly,  although  the  NH3 
mass  fractions  of  around  0.  5 are  what  one  would  expect,,  the  CONTAM  mass 
flux  is  again  overpredictod  by  5.  0 at  106°K,  the  temperature  at  which  one  expects 


the  deposit  to  be  mostly  NH^.  Since  mass  flux  is  proportional  to 
T (le.  , pV'^(T  S(T*'^s  ),  uncertainties  in  temperature  predictions  cannot 

reasonably  explain  these  large  factors.  In  section  4.4  it  was  shown  that  MULTRAN 
incompletely  expands  the  plume.  Thus,  this  routine  would  be  expected  to  concentrate 
the  mass  near  the  centerline.  Therefore  MULTRAN  is  the  most  likely  cause  of 
the  overpredictions  of  the  mass  deposition  rate. 

5.  2 Evaporation  Rate  Prediction  Errors 

As  in  the  mass  deposition  discussion,  there  are  two  possible  error  sources 
for  the  evaporation  rate  predictions:  the  equations  used  to  calculate  the  vapor 
pressures  and  the  evaporation  rate,  and  the  surface  temperature  used  in  the 
calculation.  There  is  little  question  that  the  NH3  and  N2H4  vapor  pressure  equa- 
tions are  not  extremely  accurate.  However,  the  strong  function  of  temperature 
on  vapor  pressure  makes  the  temperature  used  the  most  critical  aspect.  Like- 
wise, the  temperature  effect  all  but  overrides  any  errors  in  the  evaporation  rate 
equation  itself.  It  is  quite  likely  that,  due  to  the  kinetic  energy  impacted  to  the 
surface  by  the  impacting  gases,  the  effective  temperature  for  evaporation  would 
be  higher  than  that  recorded  by  JPL.  This  higher  temperature  would  cause 
evaporation  rates  to  be  greater  than  predicted.  If  one  examines  Figure  7 it  can 
be  observed  that  the  JPL  data  are  taken  quite  close  to  the  break  points  in  the 
curve,  those  points  where  evaporation  rate  balances  mass  flux.  Higher  tempera- 
ture of  evaporation  wovild  shift  the  JPL  data  to  the  right,  thus  providing  better 
agreement  between  experimental  and  calculated  net  deposition. 

b*  ^ Evaluation  of  CONTAM  Routines  and  Procedures 

Based  or  the  above  discussion  (5.  1),  each  routine  and  analysis  technique 
will  he  assessed. 

The  N2H4  routine  predicts  the  chamber  conditions  which  are  to  be  used  by 
MULTRAN  and  the  mass  fractions  used  by  KIN  CON.  With  the  exception  of 
hydraiiine,  the  mass  fractions  of  the  species  are  predicted  accurately.  The 
temperature  determined  by  N2H4  will  have  an  effect  on  the  results  ana  there 


were  some  problems  in  this  calculation  (see  section  3.  2).  However,  this  effect 


will  be  of  a second  order.  The  chamber  pressure  predicted  by  N2H4  agreed 
well  with  the  JPL  data.  The  N2H4  code  is  not  the  weak  link  in  the  chain.  Since 
N2H4  is  based  on  a NASA  performance  code  verified  during  development,  this 
agreement  is  to  be  expected. 

The  MULTRAN  routine  generates  the  streamlines  for  use  in  KINCON.  The 
lack  of  a proper  plume  expansion  in  MULTRAN  is  probably  a major  reason  for 
overprediction  along  the  centerline. 

The  KINCON  code  appeared  to  perform  well,  No  serious  irregularities  were 
encountered.  However,  no  kinetics  or  condensation  were  employed, and  so  KINCON 
was  not  tested  to  its  fullest. 

The  use  of  evaporation  rates  to  adjust  the  deposition  rates  worked  very  well, 
While  better  NH3  and  N2H4  vapor  pressure  data  would  be  useful,  the  order  of 
magnitude  agreement  indicates  that  the  technique  is  acceptable.  The  actual 
temperature  of  the  surface  plays  a significant  part  in  the  evaporation  rate,  and 
an  accurate  knowledge  of  this  information  is  essential  to  proper  net  deposition 
rate  calculations.  Errors  in  this  value  are  another  major  source  of  overprediction. 


6.0  CONCLUSIONS 


As  a result  of  the  comparisons  performed  on  this  effort,  the  following 
conclusions  as  to  the  effectiveness  of  the  CONTAM  code  to  predict  net  deposition 
can  be  made; 

1.  The  subroutines  of  the  CONTAM  code  (N2H4,  MULTRAN,  and  KINCON) 
coupled  with  a simple  evaporation  rate  model,  predicted  net  deposition  rates 
along  the  plume  centerline  which  were  consistently  higher  (by  factors  of  2 to  8) 
than  the  measured  rates.  Thus,  the  use  of  CONTAM  would  lead  to  ovorprediction 
of  the  centerline  deposition  rates.  The  CONTAM  predictions  did  follow  the  same 
QCM  temperature  trends  as  the  JP.L  data,  however,  and  could  be  used  to  make 
order  of  magnitude  predictions  - again  for  the  centerline  only, 

2.  MULTRAN  does  not  properly  handle  the  expansion  of  the  plume.  The 
lack  of  a proper  boundary  layer  treatment  is  a known  deficiency  in  the  code  and 
is  most  likely  responsible  for  the  overprediction  along  the  centerline,  Use  of 
the  CONTAM  code  should  give  conservative  estimates  along  the  centerline,  but 
care  should  be  taken  when  calculating  off  centerline  deposition.  Unfortunately, 
these  regions  are  by  far  the  most  critical  regions  for  contamination  considerations, 

.1,  The  N2H4  routine  appeared  to  calculate  the  mass  fractions,  with  the 
exception  of  hydra/.inu,  and  temperatures  adequately  enough  for  input  into 
MU LTRAN  and  KINCON.  The  hydrazine  prediction  would  be  important  only  in 
the  case  cf  surfaces  between  ~144°  and-'-lSO^K  or  if  the  hydrazine  flux  was 
abnormally  high,  as  with  an  aged  thruster.  While  the  lack  of  an  aging  mode  is  a 
serious  consideration  for  performance,  it  is  only  of  minor  importance  for  mc>nt 
contair.ination  assessments. 

4,  The  KINCON  model  was  not  assessed  in  its  full  options,  (condensation 
and.  vith  kinetics).  Condensation  may  play  a critical  part  of  the  hydrazine  plunie; 
kinetics  do  not.  As  used,  the  KINCON  model  appeared  to  yield  consistent  and 
reasonable  results. 

f5.  The  use  of  an  evaporation  .-ate  model  to  calculat'  net  deposition  appears 
adequate.  The  strong  effect  of  temperature  on  evaporation  points  out  the  need  to 
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know  quite  accurately  the  temperature  of  the  surface  on  which  one  is  trying  to 
predict  deposition.  This  effect  raises  some  questions  regarding  the  use  of  QCM 
temperature  to  characterize  the  evaporation  rate,  since  the  temperature  of  the 
deposit  may  be  higher  due  to  the  kinetic  energy  of  the  impinging  gases. 


7.0  FUTURE  STUDIES 


Based  on  the  results  of  the  analysis,  a great  many  questions  have  surfaced 
concerning  the  modeling  of  plume  contamination  from  hydrazine  thrusters. 

Because  the  JPL  measurements  were  made  only  every  120  seconds,  the  effects 
of  a single  pulse  and  the  off-time  could  not  be  separated.  These  measurements 
are  extremely  critical  during  the  initial  operation  of  the  thruster.  Also,  such 
measurements  would  help  determine  the  quantitative  acceptability  of  the  evapora- 
tion rate  and  the  averaging  techniques, 

To  answer  these  questions,  a test  program  is  being  conducted  at  the  Arnold 
Engineering  and  Development  Center  (AEDC).  Using  the  same  thruster  tested 
at  JPL,  AEDC  will  measure  composition  and  densities  of  the  plume  as  well  as 
net  deposition  as  a function  of  time.  Deposition  rate  measurements  will  be  made 
on  a pulse  by  pulse  basis  so  that  the  effects  of  deposition  and  evaporation  can  be 
separated.  Off-axis  and  backflow  measurements  of  the  mass  deposition  rates  for 
the  plume  expansion  are  also  being  made.  In  addition  to  this  program,  the  AFRPL 
will  be  conducting  an  analytical  program  to  improve  the  CONTAM  model,  as  identi- 
fied in  this  study,  and  to  compare  predictions  from  the  improved  code  with  the 
AEDC  test  data. 

For  contamination  predictions  alone  it  is  doubtful  that  additional  onginci 
modeling  would  be  in  order.  For  performance  predictions,  however,  there  is 
vast  room  for  improvement.  The  AFRPL  is  contemplating  a program  to  improve 
existing  codes  for  performance  and  life  prediction.  The  model  would  be  useful 
to  contamination  prediction  for  those  conditions  in  which  hydrazine  flux  is 
important.  Also,  with  a model  which  could  predict  engine  aging  (which  N2H4 
cuiTw-.ntly  can  nov.);  one  would  be  able  to  assess  whether  hydrazine  related  con- 
tamination would  be  a problem  anytime  during  the  entire  mission  life  of  the 
thruster , 
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APPENDIX  A SAMPLE  CALCULATIONS 


The  purpose  of  this  appendix  is  to  present  a sample  calculation  in  order  that 
one  may  follow  the  procedure  used  to  calculate  the  net  deposition  predicted  from 
the  CONTAM  code  coupled  with  the  simple  evaporation  model.  The  case  used  is 
that  of  the  baseline  case  (Table  A-I  summarizes  this  case);  references  are  to 
equations  or  Tables  found  in  the  text. 

TABLE  A-I 
BASELINE  CASE 


Mass  Fractions  (From  N2H4)! 

NH3  = . 51354 
N2  = .46365 

H2  = . 01571 

N2H4-  4x10-^ 

H2O  ^ .0071 

Duty  Cycle;  0.  1 sec  on/ 10  sec  off 
Distance  to  QCM;  113.  4 cm  (centerline) 

From  the  KJ.NCON  output,  density  and  velocity  as  function  of  time  can  be 
obtained.  These  values  are  used  in  equation  (1)  to  calculate  a mass  flux  per 
unit  area  as  a function  of  distance  and  then,  using  equation  (?.),  a mass  f'ux  per 
unit  solid  angle.  For  the  baseline  case  “ • 822  g/ster-sec  This  value  is 
then  converted  to  a mass  flux  per  unit  area  at  the  QCM  by  solving  equation  (2) 
for  r - 113  4 cm. 

-5  2 

m - 6 39  X 10  g/cm  - sec 

This  value  is  the  flux  predicted  to  hit  the  QCM  during  the  time  when  mass  is 
being  discharged  from  the  thruster.  The  total  buildup  rate  on  the  QCM  has  to 
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include  the  off  time  as  well.  Thus  the  value  is  adjusted  by  a ratio  of  the  "on" 
time  to  the  total  time/pulse. 

rh  = 6.  39  X 10”^  x =9.50x10“'^  g/ cm^  - sec 

Note  that  0.  1 seconds  is  the  on  time  for  the  pulse.  However,  0.  15  is  used 
becau.se  the  mass  flow  rate  at  times  between  0.  1 and  0.  15  seconds  was  con- 
sidered great  enough  to  contribute  to  the  deposition.  If  0.  1 was  used,  the 
average  value  of  the  properties  would  be  higher,  and  thus  this  effect  should  be 
minimal.  Dividing  by  10.  1 is  critical,  however,  since  the  JPL  measurements 
did  not  differentiate  between  on  and  off  time. 

To  calculate  the  individual  mass  flux  for  each  species,  equation  (3)  is  used. 
Table  A-II  presents  the  values  for  each  species  considered. 

TAB.LE  A-II 

SPECIES  DEPOSITION  RATE 

NHj  4.86  X lO""^  g/cm^  - sec 

H2O  = 6,  74  X 10"*^  g/cm^  - sec 

N2H4  = 3.  80  X 10"'^®  g/cm^  - sec 

Trie  values  in  Tabic  A-II  are  the  fluxes  that  would  be  expected  to  deposit  on 

the  crystal  if  there  were  no  evaporation.  The  following  tables  present  the  results 

of  tVis  evaporation  rate  calculations, 

TABLE  A-in 

VAPOR  PRESSURE  CALCULATIONS*  (torr) 


Species 

106 

Temperature  (K*-') 

144  172 

200 

NH3  (eq  5) 

4. 15  X 10’^ 

4.  75  X 10'^ 

-- 

H2C  (eq  6) 

2. 02  X 10'^® 

6.  35  X lO"*^ 

6.  90  X 10'^ 

8.  13  X 10 

N2ri4.  (eq  7) 

-- 

2.  37  X 10‘^'' 

1.0  X .’0’^^ 

7.  98  X 10 

* shown  only  for 

vapor  pressure  (Py)  of  interest!  10”^ 

>Pv>10-^® 

A'-IO 


TABLE  A-IV 


EVAPORATION  RATE  CALCULATIONS  » (eg  4)  g/cm^  - aec 

Species  Temperature  (°K) 

106 144  172  ~ 200 

NH3  9.70x10"®  9.51x10-4 

HgO  4.85x10"^’^  1.31x10“^°  1.  11x10"^  1.42x10"^ 

N2H4  --  --  2.74x10-^2  1.86x10-® 

* Shov/n  only  for  Evaporation  rate  (Ev)  of  interest:  10-^>Ev  >10"^'^ 

By  subtracting  the  evaporation  rates  (Table  A-IV)  from  the  deposition  rates 
(Table  A-III)  the  net  deposition  rates  are  obtained  (Table  A-V).  Note  that  a 
negative  value  for  a species  represents  no  deposition.  The  total  deposition  rates 
(sum  of  the  species)  are  then  used  to  compare  the  CONTAM  predictions  to  the 
JPL  measurements. 


TABLE  A-V 

NET  DEPOSITION  RATES  (g/cm^  - sec) 


Species 

106 

Temperature  (°K) 

144  172 

200 

NHj 

3,  90  X 10“’^ 

0 

0 

0 

H^O 

6.  74  X 10-*^ 

6. 60  X 10-^ 

0 

0 

N2H4 

3.  80  X 10-10 

3. 77  X 10-10 

3. 79  X 10-10 

0 

■TOTAL 

3,  97  X lO-"^ 

6. 98  X 10-9 

3. 77  X 10-10 

0 
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APPENDIX  B 


NET  DEPOSITION  PREDICTIONS  FOR  HYDRAZINE  THRUSTERS 

This  appendix  presents  a simple  method  for  estimating  deposition  rates 
from  hydrazine  thrusters  on  spacecraft  surfaces  without  implementation  of  the 
complex  CONTAM  code.  The  method  is  based  on  the  results  of  the  CONTAM 
evaluation,  the  JPL  test  data,  generalizations  regarding  hydrazine  thrusters, 
and  empirical  relationships  between  plume  centerline  and  backflow  mass  flux. 
Whenever  possible,  conservative  approximations  are  made, and  thus,  if  this 
analysis  indicates  that  deposition  will  occur,  a more  detailed  analysis  with  the 
CONTAM  code  may  be  called  for.  One  should  be  aware  that  these  estimates  are 
subject  to  many  of  the  same  qualifications  described  in  the  text  for  CONTAM. 

In  order  to  demonstrate  the  application  of  the  technique  discussed  here,  an 
estimate  of  deposition  from  a 5 Ibf  thruster  for  the  NATO  III  Satellite  will  be 
carried  out.  A more  detailed  analysis  of  the  deposition  from  this  thruster  has 
previously  been  done  using  the  CONTAM  Code,  and  thus  an  evaluation  of  this 
technique  is  available.  Appendix  A should  be  used  as  a guide  in  following  the 
subsequent  analysis. 

Mass  Flux  Calculations 

The  initial  step  in  analyzing  the  deposition  from  any  thruster  is  to  determine 
the  mass  flux  from  that  thruster.  While  this  mass  flux  is  a function  of  many 
variables  it  should  be  approximately  proportional  to  the  mass  flov/rate.  Since 
the  Isp  of  hydrazine  thrusters  lies  between  100  - 230  seconds,  the  flowrate  will 
var'y'  from  thruster  to  thruster  by  no  more  than  a factor  of  2.  3 times  the  thrust. 
For  the  0.  1 Ibf  engine,  the  mass  flux  ranged  from  . 7— *1.  0 g/ster-sec.  Using 


(11)  L.  P.  Davis  and  I.  L.  Witbracht,  Thruster  Contamination  Predictions  for 
NATO  III  Satellite.  AFRPL-TR-75-67,  December  1975. 
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1. 0 g/ster-sec  and  the  2.  3 factor,  the  following  can  be  used  to  calculate  the 
mass  flux! 

(Is*)'!  ^ " 23  I<T  g/sec-ster  (B-1) 

mass  flux  of  thruster  (g/sec-ster) 

F'j’  = Thrust  of  thruster  (Ibf) 

1.  0/0.  1 = ratio  of  mass  flux  of  0.  1 Ibf  thruster  to  its  thrust 

2.  3 s factor  to  account  for  low  impulse  of  the  shorter  pulses. 

Mass  Fraction  Analysis 

To  convert  total  mass  flux  to  that  for  individual  species^  the  mass  fraction 
of  the  species  in  the  plume  must  be  known.  Since  all  hydrazine  thrusters  have 
similar  performance  parameters,  the  mass  fractions  of  water,  hydrazine,  and 
ammonia  can  be  estimated  vdth  reasonable  accuracy.  For  water,  its  mass 
fraction  in  the  propellant  should  be  used;  the  highest  permissable  value  for  mono- 
propellant grade  lo  . 01.  A typical  mass  fraction  of  0.  5 for  NH3  can  be  used. 
Typical  hydrazine  mass  fractions  range  from  10“^  to  10"^,  Thif!  value  is  by  far 
the  most  difficult  to  generalize,  since  it  depends  greatly  on  duty  cycle  and  catalyst 
bed  age.  For  this  analysis,  10'^  will  be  used. 

One  could  extend  the  analysis  to  include  the  presence  of  other  contaminants, 
aniline  for  example,  but  further  analysis  would  necessitate  a knowledge  of  the 
forr."  and  ma.ss  fraction  of  the  species  in  the  plume  and,  hence,  v/ill  not  be  done 
nere  Combining  the  value  of  the  total  mass  flux  (B-1)  with  the  speciss  mass  fluxes 
g.vds  mass  deposition  rates  per  solid  angle  for  the  three  species  considered  as 
fo;lov,si 


where 


idm\ 

idfi/i  - 23  Xi  Fx  (B-2' 


whe  rc! 


Wfi/  , 

Xi  . 


= mass  flux  per  solid  angle  for  species  i 

mass  fraction  per  species  i (see  Table  B-1  for  suggested  Va.luet,) 
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TABLE  B-I 


SUGGESTED  VALUES  FOR  SPECIES  MASS  FRACTIONS 


Species 

Hydrazine 

NH3 

Water 

Other 


Suggested  Range 
10“^  - 10-4 
0.  5 

. 005  - .01 

As  determined 
by  used 


Used  in  Analysis 
. 01 
0.  5 
0.  01 

Not  considered 


Estimation  of  Species  Flux  to  Spacecraft  Surface 

There  are  four  factors  in  converting  the  flux  per  solid  angle  (eq.  B-2)  to 
actual  fluxes  at  the  spacecraft  surface:  distance  from  the  exit  plane  to  the 
surface  (r),  the  angle  between  the  thruster  centerline  and  a line  drawn  from  the 
exit  plane  to  the  surface  (0),  the  inclination  of  this  surface  with  respect  to  that 
line  (^),  and  the  duty  cycle  to  be  considered.  The  geometric  parameters  arc 
depicted  in  figure  B-I. 


(cirTi\ 

dir)  i 

each  species  i by  r2.  This  converts  the  units  to  g/cm^  - sec  at  that  distance. 

The  angular  effect  (0)  cannot  be  calculated  well,  especially  for  the  backflow 
region.  Thus,  the  flux  at  any  angle  6 is  estimated  as  a fraction  of  the  centerline 
flux.  Figure  B-2  gives  the  function  used  to  estimate  this  fraction.  For  relatively 
small  value, 8 of  0,  the  theoretical  Hill  and  Draper  approximation  is  used.  For 
larger  angles,  the  set  of  measurements  performed  at  JPL  on  an  expansion  of 
nitrogen  gas  through  a nozzle  is  used.  These  measurements  are  now  felt  to 

be  h.'.gh  m the  backflow  region^  but  their  use  represents  a conservative  estimate. 


\12)  J.  E.  Chirlvella,  Molecular  Flux  Measurements  in  the  Backflow  Region  of 
a Nozzle  Plume.  NASA  Tech  Memo  33-620,  July  1973 

(13)  R.  Passa.maneck,  Personal  Communication 
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THRUSTER/VALUE 

ASSEMBLY 


Figure  B-1.  Geometrical  Baaia  for  Backflow  Equations 
(From  Reference  1.2) 
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f(tf)  = (dm/dfil  /(dm/dfl) 


I'  ’ 

I ' 


I Figure  .0-2,  Measured  Values  of  {($)  vs 

. (From  Reference  12) 


<1 

ii 


EMPIRICAL  CURVE 
BASED  ON  BACKFLOW 
MEASUREMENTS 


Backflow  Angle  0 


MMLlliftiiiil. 


These  were  the  same  measurements  used  in  the  NATO  III  analysis.  The  surface 
inclination  angle  0 is  taken  into  account  with  a cos  factor  in  the  equation.  The 
flux  calculated  so  far  is  for  the  on-time  of  the  thruster,  which  should  be  reduced 
by  the  percent  of  the  duty  cycle  i.  e.  on-time/total  pulse  time. 

These  four  factors  result  in  the  following  equation  for  the  flux  through  the 
surface: 

23  Xi  Ft  f(6')  (DC)  (B-3} 

m^  a ' 

r^  cos  'p 

where:  rhi  is  the  flux  through  the  surface  in  g/cm^  - sec  per  specie, s i 
f(e)  = the  emperical  function  of  <p  from  Figure  B-2 

0 a the  angle  between  the  surface  normal  and  a line  joining  the  surface 
with  the  exit  plan(5 

r = the  dih'tance  between  the  exit  plane  and  the  surface  in  cm 
DC  = percent  duty  cycle,  on-time/total  pulse  (dimensionless) 

(DC  a 1.  0 for  steady  state) 

The  fluxes  calculated  from  equation  B-3  represent  the  saxne  fluxes  as  those 
of  Table  A-2,  i.  e.  , the  deposition  rate  of  the  species  onto  the  surface,  The  n_^ 
deposition  rate  at  the  surface  is  then  calculated  using  the  procedure's  described 
in  Appendix  A.  This  value  is  the  average  deposition  rate  during  a pulse  (both  on 
and  off)  and  s.hould  b<3  multiplied  by  the  time  of  pulsing  to  yield  total  net  deposition. 
Since  evaporation  is  taking  place  continually,  that  rate  may  then  be  applied  to  the 
off  time  to  determine  the  length  of  time  for  which  there  will  bi3  deposits  on  the 
surface. 

Sample  Case 

The  procedures  described  above  will  be  used  to  assess  the  net  deposition 
rate  from  a Ibf  thruster  used  on  the  NATO  III  Satellite.  The  parameters  to  be 
used  are  described  on  Table  B-II. 
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TABLE  B-Il 

PARAMETERS  FOR  SAMPLE  CASE 


Parameter 

Percent  Duty  Cycle 
(on/total) 

Total  Pulse  Time 

Surface  Angle 

Distance  to  surface 

Thrust 

Angle  from  centerline 


Nomenclatur  e 
DC 

r 

F 

e 


Value 
. 150 

(.  09  on/.  511) 
. 601 
60° 

30  cm 
5.  0 Ibf 
120° 


From  figure  B-2,  f(120°)  « 2 X 10"^.  Applying  equation  B-3  to  this  case 
using  the  mass  fractions  from  table  B-I  yields  the  species  deposition  rates  listed 
in  Table  B-lIl.  That  table  also  compares  the  results  with  those  obtained  for  the 
NATO  III  thruster  using  CONTAM. 

TABLE  B-III 


Species 

NH3 

H2O 

N2H4 


SPECIES  DEPOSITION  RATES  (g/cm^  - sec) 


Estimation  Technique 

CONTAM 

3.  82  X 10"® 

1.49  X 10"^ 

7,  65  X 10"’^ 

3.  87  X lO"”^ 

7.65  X 10“'^ 

9, 38  X 10"® 

It  can  be  seen  from  table  B-lII  that  the  technique  predicts  deposition  rates  which 
are  within  an  order  of  magnitude  of  values  predicted  by  CONTAM  and  conserva- 
tive. Thus  this  technique  should  give  good  first  estimate  predictions.  The  next 
step  is  to  use  the  evaporation  rates  from  the  surface  to  calculate  the  deposition 
expected  as  a result  of  a single  pulse.  Surface  temperatures  216°K  (-72°F)  and 
144°K  (••200°F)  will  be  used.  The  vapor  pressure  and  evaporation  rate  equations 
diacusBod  in  Appendix  A and  the  text  were  used  to  produce  Table  B-IV.  By 
subtracting  this  rate  from  the  deposition  rate  (Table  B-III),  estimated  net  depo- 


sition is  obtained  (Table  B-V). 
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TABLE  B-IV 


Species 

NH3 

H2O 

N2H4 


SPECIES  EVAPORATION  RATES  (g/cm-sec) 


216°K 

144°K 

461 

9. 51  X 10'^ 

2.  00  X 10“^ 

1.  31  X 10-^° 

4.  73  X 10-'7 

6.  5 X 10“®9 

Species 

NH3 

H2O 

N2H4 

Total 


TABLE  B-V 

NET  DEPOSITION  RATES  (g/cm^ 

216°K 

0 

0 

2.  92  X 10-'!' 


sec) 


144°K 


7.  59  X 
7.  65  X lO"''' 
1.  52  X 10“^ 


2.  92  X 10-'!' 

This  analysis  indicates  that  the  net  deposition  during  each  second  is  2,92  X g/cm^ 
and  1.  52  X 10"^  g/cm^  at  216°K  and  144°K  respectively.  If  the  pulse  train  is  500 
pulses  long,  then  this  deposition  occurs  for  0.601  X 500  or  300,  5 seconds  for  a 
net  mass  buildup  of  8.  77  X 10'®  g and  4.  57  X 10’*^  g at  2l6°K  and  144°K  respectively. 
Using  the  appropriate  evaporation  rate  will  permit  the  determination  of  the  length 
of  time  required  to  evaporate  the  deposition. 


Graphical  Calculations 

A graphical  method  for  estimating  the  net  deposition  can  be  performed  using 
the  following  steps  (sample  calculation  is  shown  on  graphs). 

1.  Calculate  the  quantity  [DCFq'/r^  cos^]  (for  the  sample  above  this  value 
is  1,  66  X 10-3). 

2.  On  Figure  B-3,  intersect  a horizontal  line  from  quantity  (from  1)  with  a 
vertical  line  at  the  angle  corresponding  to  the  degrees  off  centerline  (120°  for  the 
centerline).  The  mass  flux  is  determined  by  the  point  of  intersection.  (In  the 
sample  case  it  is  1.  1 x 10"^  g/cm^  - sec). 
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3.  The  mass  £lux  should  be  multiplied  by  the  mass  fraction  of  each  species 
considered.  Figure  .B-4  contains  suggested  values. 

4.  The  flux  of  each  specie  should  be  tepresented  by  a horizontal  line  on 
figure  B-4  at  the  appropriate  ordinate  value. 

5.  The  solid  lines  tabled  NH3,  H2O,  and  N2H4  on  Figure  B-4  represent 
the  evaporation  rate  of  the  species.  Thus  the  difference  between  the  intersection 
of  a vertical  line  through  a temperature  with  the  horizontal  line  (from  4)  and  the 
solid  line  (evaporation  rate)  is  the  net  deposition  for  that  species.  If  this  value  is 
negative,  there  will  be  no  deposition  of  that  species.  The  temperature  s.t  which  the 
horizontal  line  intersects  the  solid  line  is  the  temperature  above  which  there  will 
be  no  deposition  of  that  species. 

6.  For  the  example  at  2l6°K,  only  N2H4  is  predicted  to  deposit  v/ith  a rate 
of  10“^  - 3 X 10"^  = 7 X lO"”^  (compared  with  3 X 10“^  fci-  a rigorous  calculation). 
For  144°K  all  the  hydrazine  and  water  will  deposit  for  a total  of  2 X 10"^  (com- 
pared with  1.5  X 10'^  for  the  rigorous  calculation),  Ono  can  also  deduce  that  no 
NK3  will  deposit  until  128°K. 
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